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ABSTRACT 

4U 1822—371 is one of the proto-type accretion disk coronal sources with an orbital period of 
about 5.6 hours. The binary is viewed almost edge-on at a high inclination angle of 83 degrees, which 
makes it a unique candidate to study binary orbital and accretion disk dynamics in high powered X-ray 
sources. We observed the X-ray source in 4U 1822—371 with the Chandra High Energy Transmission 
Grating Spectrometer (HETGS) for almost nine binary orbits. X-ray eclipse times provide an update 
of the orbital ephemeris. We find that our result follows the quadratic function implied by previous 
observations; however, it suggests a flatter trend. Detailed line dynamics also confirm a previous 
suggestion that the observed photo-ionized line emission originates from a confined region in the 
outer edge of the accretion disk near the hot spot. Line properties allow us to impose limits on 
the size of accretion disk, the central corona, and the emission region. The photo-ionized plasma is 
consistent with ionization parameters of log £ > 2, and when combined with disk size and reasonable 
assumptions for the plasma density, this suggests illuminating disk luminosities which are over an 
order of magnitude higher than what is actually observed. That is, we do not directly observe the 
central emitting X-ray source. The spectral continua are best fit by a flat power law with a high 
energy cut-off and partial covering absorption (TVh ranging from 5.4-6.3 x 10 22 cm" 2 ) with a covering 
fraction of about 50%. We discuss some implications of our findings with respect to the photo-ionized 
line emission for the basic properties of the X-ray source. 

Subject headings: accretion, accretion disks — binaries: eclipsing — stars: individual (4U 1822—371) 
— X-rays: binaries 



1. INTRODUCTION 

Eclipsing X-ray binaries are still fairly rare as they re- 
quire a view of the accretion disk that is within sev- 
eral degrees of the edge. The low-mass X-ray binary 
(LMXB) 4U 1822-371 is a specifically rare object be- 
cause it also has a relatively short orbital period which 
produces eclipses every 5.57 h and a light curve which 
is ve ry sensitive to the st ructure of the accretion disk 
rim ([White fe Holt I Il982| ) . Repetitive features in the 
light curve as well as the fact that the eclipse is only 
partial provide some direct geometrical constraints such 
as a viewing angle i within 75° and 85° as well as a disk 
radiu s of (5 — 7.3) x 10 10 cm for this range of viewing an - 
gles (IMason fc Cordova I IT982I : IHellier fc Mason I Il989h . 
iHeinz fe Nowak I (J2001I ) determined the most accurate 
inclination of i = 82.5 ± 1.5 from Rossi X-Ray Timing 
Explorer (RXTE) data. 

Its generally hard X-ray spectrum is an indicator that 
there is an accretion disk corona (ADC) in the line of 
sight to the central X-ray source which also partially 
obsc ures the source with a radial extent of ~ 3 x 10 10 
cm (IWhite et al. I [19811 IWhite fc HolTI fl98l . X-ray 
spectra from the pre- Chandra era were complex to model 
as they seem to require additional and unresolved soft 
emission line structures in addition to a significant Fe K 
fluorescence line. This made 4U 1822 —371 the archetypi- 
cal example of a LMXB AD C source (jWhite et al. 1119811 
119971 IParmaret aT~l 120001 ) . The first Chandra observa- 
tion in 2000 finally resolved these line structures into 
discrete emission lines from photo-ion ized O, Ne, Mg, 
Si, S, and Fe ions (|Cottam et al. ||2001| ). Phase-resolved 
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spectra also suggested that the line emission originates 
from a highly confined area, likely the illuminated hot 
spot. From the persistence o f the Fe K fluor e scenc e line 
throughout the binary orbit, iCottam et al. I ()2001l ) also 
concluded that its emission must come from a more ex- 
tended region around the disk. A recent survey of all X- 
ray binaries in the Chandra archive indeed shows that in 
nearly all cases narrow line fluorescence originates from 
more or less spherically distribu ted material around the 
centrally illuminating source ()Torreion. et al. |[2010D . 

More detailed dynamical information in LMXBs are 
hard to come by simply because a clear identification 
of the low-mass com panion at distances of several kpc 
is usually impossible. Uonker fc van deFLv lis (200l|) de- 
tected 0.59 s pulsations from the neutron star in 4U 
1822—371 and determined a donor star mass of 0.4 
Mq assuming a neutron star mass of 1.4 M©. From 
pulsation time delays they also determined that the 
orbit must be circular. By observing Bowen fluores- 
cence from the X-ray heated face of the companion star, 
ICasares et al. I ( 2003) determined a radial velocity semi- 
amplitude of 300 km s , which yielded a comp anion 
mass of 0.36 M Q . However iCowlev et al. I (|2003l ) used 
strong He I absorption to measure a lower limit of the 
orbital velocity of the donor star to 234 km s" 1 , consis - 
tent with previous measurements (jHarlaftis et al. 1 11997) . 
providing a lower limit to the donor m ass of again 0.4 
M Q . In addition, ICowlev et al. I (|2003[ ) also measured 
He II and H a emission line velocities and determined a 
systemic velocity of -103 km s _1 . They argued that such 
a large negative value in direction of the Galactic center 
might indicate that 4U 1822—371 is either part of the 
halo population or it was given a significant kick when 
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the neutron star was formed. 

In this paper we use phase- and highly spectrally re- 
solved X-ray spectra over almost nine binary orbits in 
order to further investigate the origin of the recombina- 
tion line emission and use their dynamics to diagnose 
properties of the accretion disk. 

2. DATA REDUCTION 

We observed the source 4U 1822-371 on 2008 May 
20 (ID9076) and 23 (ID 9858) with the Chandra HETGS 
(|Canizares et al. I [2000) in a standard configuration via 
the HETG Guaranteed Time Observation program (see 
table 1). With a total exposure time approximately 
150ksec, these two observations cover seven (three and 
four, respectively) full binary periods. We reprocess all 
the observations using CIAO Version 4.2 with the most 
recent CALDB products. 

For each observation, we redetermine the positions of 
the zeroth-order ccntroid using the f indzo script for op- 
timal wavelength accuracy. We extract all first-order 
spectra in HEG and MEG and co-add them into one 
single spectrum. This has the consequence that we 
adopt a spectral resolution near 0.021mA provided by 
the medium energy transmission grating grid throughout 
the entire wavelength band. The sacrificed resolution is a 
trade-off for easier handling of the spectra with respect to 
signal-to-noise ratios, specifically for the phase-resolved 
spectral analysis. Light curves and s pectral analysis are 
based on ISIS Versions 1.5.0-1.6.1 (JHouck &; Denicolal 
[200(1 . 

3. ANALYSIS AND RESULTS 

Eclipses are rare in LMXBs. They, however, contain 
critical dynamical information with respect to the geom- 
etry of the disk and its environment. A broad eclipse, 
for example, would indicate an extended companion. A 
partial eclipse may indicate the presence of large scale 
scattering structures such as ADC. In a first step in the 
analysis we derive an improved ephemeris based on all 
three HETGS observations and test if the orbital period 
is changing. A second step investigates orbital phase 
variations of bright X-ray lines. In a third step we fit a 
phase-binned spectrum at an orbital phase which allows 
a direct view of the photo-ionized region with appropri- 
ate continuum and photo-ionization models. 

3.1. An Ephemeris Update 

All existing HETGS observations of 4U 1822-371 yield 
a total of 9 complete eclipse light curves. The year 
2000 data indicated a slowly increasing orbital period. 
Here we first folded each of the observat ions into phase 
bins u sing a fixed period as determined bv lParmar et alH 
(2000). Then we fitted each phase- folded light curve with 
Gaussian and sinusoidal functions. The results are shown 
in the right upper panel of Figure [1] for ID9858 as an il- 
lustration. The new eclipse time can be well-determined. 
Once we join this new ecli pse time to the alrea dy pub- 
lished X-ray eclipse times (|Parmar et al. 2000, and the 
reference therein), we can fit the timing residuals with 
respect to the best-fit linear ephemeris (x 2 = 35.99 for 
19 dof ) using the quadratic function (shown in left lower 
panel of Figure [lj. The new curve is somewhat flat- 
ter than the old one and leads to an updated quadratic 



ephemeris given by: 

Tdc = 2445615.30927(41) + 0.232109006(67) N (1) 
-0.96(18) x 10~ n V 2 

Where N is the cycle number and linear ephemeris, and 
errors are at 90% confidence on the last 2 digits. For com- 
parison, the dotted line show s the quadratic fit ob tained 
for the old X-ray datasets bv lParmar et al. I J2000). The 
new X-ray eclipse times derived from the HETG obser- 
vations are given in table [2] with uncertainties given at 
90% confidence. 



3.2. Phase Resolved Line Properties 

Observed spectral features are highly phase-dependent. 
The new observations cover almo st four times th e num - 
ber of binary orbits than what iCottam et al. I (J2001I ) 
had available, which allow us to investigate line varia- 
tions with orbital phase in much more detail. The lower 
right panel in Figure [T] shows the phase separation for 
our stronger line detections. We center the first phase 
bin at eclipse minimum. To study the H-like and He- 
like lines, except the narrow phase bin (~ 0.08) during 
eclipse, we choose a phase bin width of 0.16. This secures 
enough counts for the line features in each bin, and the 
Ne x, Ne ix, Mg XII, Mg xi, O vin, and O vn emission 
lines are clearly detected with this choice. We divide 
the orbit into 40 even, but overlapping, phase segments. 
The Fe xxvi and Fe Ka have very weak signal with this 
choice of (correlated) phase bin, consequently, we addi- 
tionally chose a set of five, non-overlapping phase bins 
as marked in upper right panel of Figure [T] They are 
denoted as upper(l), decreasing(2), bottom(3), eclips- 
ing(4), and rising(5) respectively. 

3.2.1. Fe lines 

The crude five phase binning of the Fe xxvi and Fe Ka 
lines smears out most orbital information. Detections 
are shown in Table [5j and no line shifts are found. How- 
ever, their flux trends follow the shape of the light curves, 
reaching a minimum during eclipse and a maximum at 
the light curve's maximum. For all five phases, the cen- 
troid of the Fe K line is consistent with 1.93 A which 
covers wavelengths of Fe I to Fe XX from the cold, near- 
neutra l medium, similar to the finding by ICottam et alTI 
(2001). However, we find a lower limit to the line broad- 
ening of about 370 km s _1 . 

3.2.2. Bright Line Dynamics 

The preliminary analysis in ICottam et al. 1 1)20011 ) only 
had two binary orbits at hand, which provided four phase 
bins and tentative suggestions with respect to an orbital 
phase dependence of bright line centroids. Incorporat- 
ing the seven orbits in the new observations allows for 
a much more detailed analysis of the line phase depen- 
dence. These line properties include line centroids, line 
widths, and line fluxes. The results are shown in Fig- 
ures [5] (for Lya lines) and [6] (for intercombination lines) 
and Tables 3 and 4. These are the brightest lines in the 
X-ray spectrum. There are other weaker lines present in 
the spectrum; however, we perform the phase dependent 
analysis only for the bright lines. 
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All of the lines show significant wavelength shifts with 
orbital phase. Most lines follow a very similar shift pat- 
tern. At phase the lines appear highly red-shifted and 
exhibit values between 400-500km s~ . These values de- 
cline to zero near about phase 0.25 and turn into increas- 
ing blue-shifts of up to 400-500km s _1 near phase 0.5 
where in most cases the shifts reverse again. The most 
prominent lines showing this pattern are Ne ix, Mg XI, 
Mg xii, and Si xiv. In Ne Xthe shift reversal happens 
later at phase 0.8, at O vnand O vin somewhat earlier 
than mid-phase. The wave pattern appears most clearly 
in Ne ix, Mg XII, and Si xiv, which are also the brightest 
in the sample. Uncertainties vary between 100 and 250 
km s" 1 depending on the phase and the detected flux. 
This pattern is very consistent with a spatially confined 
emission region moving with orbital phase. We illustrate 
this with a cartoon structure of the system and the phase 
light curve highlighting the phase bin location as well as 
the line spectrum for the case of Ne ix in Figure [71 
Shown phases include 0, 0.1, 0.42, 0.68, 0.83, and 0.92. 

The second part in Figure [5] also shows the veloc- 
ity width of the lines. Even though there appears to 
be some subtle variations with phase, in most cases 
the widths appear rather stable over the entire or- 
bit. The velocity widths for all lines appear very sim- 
ilar near 400 km s _1 with uncertainties of around 150 
km s~ per phase bin, again depending on detected flux. 
The very similar widths are consistent with a radius of 
(1.14±0.7) x 10 n cm. 

Table [7J and Figure |4] list the results of our detections 
for He-like lines of O VII, Ne ix, and Mg XI. The inter- 
combination line flux is very bright and appears broad, 
while the resonance line is weak, and we do not detect 
the forbidden line component. The resulting G ratios, 
defined as (i+f)/r, are all consistent with a pure recom- 
bination plasma. The R ratios, defined as f/i, are all very 
small indicating either high levels of photo-excitation or 
the presence of a high density plasma. To remove the or- 
bital smear we compute the R and G values for the five 
phase bin case, which are shown in Figure [3] and Table El 
In orbital phases 2 and 4 we do not have enough statis- 
tics to compute values, however for phases 1, 3, and 5 
we compute values and value limits which are consistent 
with each other. 

3.3. The Photo-ionized X-ray Spectrum 

We fit spectra for the five phase bins described in Sect. 
3.2 and find that it can be fit well at most phases us- 
ing a cut-off power law plus a soft black body compo- 
nent. The norm variations in both components follow 
the flux change in the light curve. The power law in- 
dex is very flat with values in the range w —0.4 — 0.2. 
This flat power law reflects the fact that we likely observe 
indirect emission affected by a central hot corona as al- 
ready found in RXTE and ASCA data (jHeinz k. Nowak I 
2001). The blackbody component has an average tem- 
perature of 181±11 eV with slight variations, which pos- 
sibly reflects some correlation with the photo-electric ab- 
sorption function in the fit. We detect some significant 
absorption using the updated Tbnew function in Xspec 
of (1.6±0.3) x 10 21 cm~ 2 throughout all phases, which 
is nea r the upper limit presented by iHeinz fc Nowakl 
(|2001f h The highest blackbody flux in the light curve 



amounts to 3xl0 34 erg s _1 , which we estimate is only 
about 2% of the total flux (see below) and provides an 
emission radius of about 20 km. Figure [8] (left) shows 
the continuum fit to phase 1 with blackbody, power-law, 
and gaussian line components. 

Alternatively i t ha s also been suggested by 
IHeinz fc Nowak I (|2001f ) and others that the spec- 
tra should be fit using the power law with additional 
partial covering absorption. We have simultaneously 
fit phases 1-5 with such a function wherein we forced 
the power law slope to be the same for all phases (the 
fit yields T = 0.656 ± 0.002, with 90% confidence error 
bars), but let the power law normalization, line param- 
eters, partial covering fraction, and covering column 
be free parameters for each phase. Figure [8] (right) 
shows the results of this fit for the phase 1 spectra, 
with the photo-ionized lines represented by gaussian 
functions. We find that this model simultaneously fits 
all 5 spectra very well (x 2 = 2266.6/2176 DoF).' The 
direct and covered power law represent a normalized 
flux that ranges from 0.57-1. That is, at least 43% 
of the total power law flux is completely blocked in 
the faintest phase, i.e., phase 4. (In all likelihood, a 
larger fraction of the ADC flux is completely blocked, 
as phase 1 is unlikely to be completely uncovered by 
the disk rim.) The partial covering fraction varies 
between 49-52%, while the partial covering column 
varies between 5.4-6.3 x 10 22 cm -2 . In this case the ISM 
absorption remains very low at (2.3±0.3) x 10 20 cm~ 2 . 

The line emitting region is best viewed during phases 
1 and 5, with the least dynamical smear during phase 5. 
During phase 1 the line emitting region is more directly 
viewed, but here the lines switch from red-shift to blue- 
shift and the region is likely already partially obscured 
by the central corona. Figure [3] thus shows the Ne ix line 
strongest and most narrow. We select phase bin 5 for the 
more detailed spectral modeling. The main goal of this 
analysis is to get an estimate for the ionization param- 
eter of the emitting region. Photo-ionization properties 
are generally described by the parameter £ = L/(n*r 2 ), 
where L is the source X-ray luminosity, n is the electron 
density, and r is the distance to th e illuminating X-ray 
source ()Kallman k. McCrav 1119821 ). From Figure [5] we 
already see that the O VII and O vm line strengths are 
similar, which is difficult to achieve with a single ioniza- 
tion parameter. We thus expect a range of parameters. 

We want to keep this analysis as simple as possible 
and use the new embedded Xstar function photemis in 
Xspec. Fitting the lines with a single ionization param- 
eter, we obtain log £ = 2.6 [erg cm/s], which fits the 
H-like ion lines very well but not the He-like ion species. 
However, even though a second ionization parameter of 
log £ = 1.9 [erg cm/s] would fit these lines, the fit itself 
is not acceptable as it would produce significant radia- 
tive recombination continua (RRC) for O vm, Ne IX and 
Ne x, which we do not observe. A viable solution to make 
up for this deficiency is to assume a shallow tail of lower 
ionization parameters (i.e. below log £ = 2.6 [erg cm/s]), 
allowing for high enough plasma temperatures to smear 
out the RRCs but produce enough line flux to fit the He- 
like lines. The range of this tail depends on abundance 
choices and here we cannot find unique solutions. 

4. DISCUSSION 
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4U 1822—371 is known as one of the prototypical ADC 
systems and shows a strong orbital dependence of its X- 
ray properties. Our coverage of many orbital periods 
allows us to study these dependencies in detail. The 
updated eclipse times from the three Chandra observa- 
tions are not inconsistent with the tre nd provided by 
previ ously published X-ray eclipse times (jParmar et al. I 
2000). But even though our result does follow the 
quadratic function implied by these previous observa- 
tions, it does suggest a flatter trend. This is still con- 
sistent with a previous result, which provided the last 
updated improved ephemeris of 4U 1822—371 d erived 
from UV/optical light curves (|Bavless et al. 1 120091) . But 
in their analysis, the quadratic term in the optical is 
also c onsistent with the one in X-ray by iParmar et al. I 

(2000) within the measurement errors, and the eclipse 
time in the o ptical lags behind by about 100 seconds. 
More recently, IBurderi et al. I (|2010t) revisited X-ray ob- 
servations o ver the last 30 y rs co nfirming the quadratic 
function by IParmar et al. I (|2000f ) , but did not include 
recent optical and UV data. In their analysis the eclipse 
times from our Chandra observations showed much larger 
delays which appear closer to the quadratic function by 
IParmar et al. I (|2000l ) . however the delays also remained 
significantly below the function. One possible explana- 
tion for the discrepancy with our result is that they opted 
not to fold the light curves, which we think is warranted 
given the multiple eclipses within these data, and the fact 
that we see intermittent dips in the spectra, some occur- 
ing quite close to the main eclipses. In any case, the trend 
solidifies, that the orbital period change in this system is 
larger than expected from simple magnetic braking and 
losses due to gravitational radiation. 

The line variations with orbital phase can probe geo- 
metrical properties such as accretion disk size, ADC size, 
and also a phase-resolved ionization balance of the line 
emitting region. While most of the lower Z ions show sig- 
nificant dynamical properties along the binary orbit with 
respect to all of their line properties, the iron lines only 
show flux changes consistent with the changes observed 
in the light curves. The Fe Ka fluorescence line is visible 
at all phases. Its flux is strongest opposite of the eclipse, 
when the disk is fully exposed to the observer. During 
eclipse the fact that the line flux is still about 30% of its 
maximum indicates that the line emitting region extends 
significantly beyond the size of the companion. We also 
do not believe that the illuminated surface of the stellar 
companion is a major source of Fe fluorescence, because 
at an edge-on view and opposite to eclipse the disk would 
likely block the flux and we would expect maxima during 
the decreasing and bottom phases (Table Et - Similarly , 
using fluorescence probabilities calculated by iBai I (|1979l ) 
we also expect maxima at these phases and not at the op- 
posite (upper) pha se. Instead we a rgue, consistent with 
the suggestions by iCottam et al. I (|2001|) . that the fluo- 
rescence materials c ome from an ext e nded region above 
the disk. In fact, as lTorreion. et al~l (|2010D showed in a 
large survey of Fe K fluorescence in X-ray binaries includ- 
ing 4U 1822—371 that the emissions come from a large, 
spherically distributed volume around the accretion disk. 
The width of Fe Ka fluorescence line are also better con- 
strain ed than the previous detection by ICottam et al. I 

(2001) and we determine a lower limit of the emission 



radius of R > 1 x 10 95 cm. We do not know how well 
this fi ts into the previous conclusion by iBavless et al. I 
(2009) with respect to the existence of a layer of cooler 
material at the base of the wind seen in the UV, but the 
X-ray data do not exclude this possibility. 

We also detect weak line emissions from hot 
Fe xxvi ions. As is observed for the cool Fe Ka fluo- 
rescence, the Fe xxvi, Lyman a line does not show any 
line shifts and broadening with respect to orbital phases. 
Its flux roughly follows that of the light curve and is be- 
low detection during eclipse. This indicates that the line 
originates from a central region within the disk radius, 
likely from parts of the central corona but not the disk it- 
self. Its emission volume is then limited by the size of the 
companion duri ng eclipse. Its Roche- lobe filling radius is 
about 0.54 i? Q (jBavless et al. Il2009f ) and thus the emis- 
sion size is less than 3.5 xlO 10 cm and thus about 0.65 
tdisk assuming an orbital separation of a = 1.33 xlO 11 cm 
and a disk radius 0.4 a (values taken from IBavless et aLl 
l2009h . 

Photo-ionized X-ray line e mission in 4U 1822 — 371 
has been observed previously (|Cottam et al. ||2001D and 
phase-resolved spectra suggested that the emissions from 
line recombination originate in an X-ray illuminated 
bulge located at or near the predicted point of impact 
of the accretion stream from the companion and the 
disk. Our observations over several binary orbits not 
only confirm this prediction, but allow us to study this 
phenomenon in detail. The lower Z line emissions clearly 
show red- and blue-shifted motion with orbital phase 
which can only be interpreted by motion of a localized 
line emitting area at the hot spot. The maximum abso- 
lute velocities from the line shifts of the most prominent 
lines (except Fe xxvi, see above), i.e., O VII, O vm, 
Ne ix, Ne x, Mg XI, Mg XII, Si xiv, yield a velocity of 
500 km s _1 , which is consistent with the orbital speed 
projection expected from the orbital parameters of 4U 
1822—371 at the inclination of 83°. Using 40 overlapping 
phase bins we observe the lines moving at all times and 
variations in the line widths appear very small, which 
restricts the size of the emitting area to < 10 6 cm. In 
fact, the very similar width in all lines places the emit- 
ting region at a consistent distance of 1.1 x 10 n cm to 
the source, which is about twice the actual outer disk 
radius and more close to the separation of the two stars. 
However, even though the uncertainties in the width in- 
clude the outer disk radius we predict that there is local 
turbulence involved in the width. 

There are a few more interesting implications from 
the broadband analysis of the X-ray spectrum. The 
fact that we fit a soft blackbody spectrum to the soft 
part of the X-ray spectrum is consiste nt with one solu- 
tion o btained from the ASCA analysis (|Heinz fe Nowakl 
|2001[ ): however, the blackbody temperature we measure 
is very much lower. The properties we determine pro- 
vide an emission radius close to 20 km, which appears 
like the emission from the neutron star itself. However, 
the line e mission properties imply a radically different 
picture. As lHeinz fe Nowakl ([20011 ) already suspected in 
their analysis of the ASCA data, we likely do not observe 
the X-ray source directly and the true source flux should 
be considerably higher. The ionization parameter of the 
line emitting region is directly related to the illuminat- 
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Table 1 

Observations for 4U1822-371 






obsID 


MJD interval 


Obs. Start 


Exp.(ks) 


PI 


671 

9076 

9858 


51779.681 - 51780.163 
54606.949 - 54607.708 
54609.551 - 54610.524 


2000-08-23 16:20:37 UT 
2008-05-20 22:46:21 UT 
2008-05-23 13:14:05 UT 


39.95 
63.68 
82.19 


Kahn, M. 
Canizares, C. 
Canizares, C. 



ing source flux (see Sect. 13. 3[) . At values between log £ 
= 1.9—2.6, a disk radius of 5xl0 10 cm, and a plasma 
density between 10 12 and 10 14 cm -3 , for the thin illumi- 
nated layer the illuminating luminosity has to be much 
larger than 10 37 erg s _1 . This also has a consequence 
for the blackbody emitting radius, which once adjusted 
to the higher source luminosity is more of the order of 
> 100 km. This makes the solution with the blackbody 
component highly unattractive as it removes the possi- 
bility that we see the neutron star, but it is also rather 
weak to account for inner disk emissions. In this respect 
we argue for the solution with the partial covering com- 
ponent, i.e., the X-ray source is a power law in nature, 
likely scattered into our line of sight over the spatially 
very extended corona, and due to the edge-on view we 
observe some of this emission both blocked and heavily 
absorbed by the disk rim. 



• Line widths appear larger than the outer disk ra- 
dius and likely involve local turbulence. 

• The line emitting region exhibits high ionization 
parameters, while missing RRCs constrain values 
to be above log £ > 2. This strongly suggests that 
we do not directly observe the X-ray source. 

• The luminosity implied by the range of ionization 
parameters favors an emission model which consists 
of a cut-off power law plus a 50% partial covering 
absorption. 



We thank Mike Noble, John Houck for their help with 
ISIS. We gratefully acknowledge the financial support of 
Chandra X-Ray Observatory archive grant AR0-11005X. 



5. CONCLUSIONS 

The orbital phase-resolved analysis of the X-ray photo- 
ionized region in 4U 1822—37 obtained the following re- 
sults: 



• The orbital X-ray ephemeris update suggests a flat- 
ter trend with respect to the previously proposed 
quadratic change function. 

• Fe xxvi line emission arises from a central corona 
of size < 3.5 x 10 10 cm. 

• Lower Z line emissions show orbital shifts consis- 
tent with a small local region at the rim of the disk, 
exhibiting orbital blue- and red-shifts of the order 
of 500 km s" 1 . 

• The line emitting region has to be the inner illu- 
minated side of a bulge at the hot spot because 
it becomes entirely eclipsed when the outer side of 
the spot faces the observer. 

• The latter fact also implies that the illuminated 
layer on top of the disk is thin, at least thin enough 
not to exceed the outer disk rim. 
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Table 2 

New X-ray eclipse times for 4U1822— 371 with 
uncertainties given at 90% confidence. 



obsID 



JDfT 



© 



Uncertainty Cycle 



671 

9076 

9858 



2451780.36424 
2454607.69102 
2454610.24377 



0.00079 
0.00045 
0.00037 



26561 
38742 
38753 



X-Ray Lines in 4U1822-371 



Table 3 

Lya lines along the phases for 4U1822— 371 a 



phase 



A 



SiXIV Lya 
Flux 



A 



MgXII Lyo 
Flux 



NcX Lya 
Flux 



OVIII Lya 
Flux 



1 0.0 

2 0.025 

3 0.05 

4 0.075 

5 0.1 

6 0.125 

7 0.15 

8 0.175 

9 0.2 

10 0.225 

11 0.25 

12 0.275 

13 0.3 

14 0.325 

15 0.35 

16 0.375 

17 0.4 

18 0.425 

19 0.45 

20 0.475 

21 0.5 

22 0.525 

23 0.55 

24 0.575 

25 0.6 

26 0.625 

27 0.65 

28 0.675 

29 0.7 

30 0.725 

31 0.75 

32 0.775 

33 0.8 

34 0.825 

35 0.85 

36 0.875 

37 0.9 

38 0.925 

39 0.95 

40 0.975 



6 - 185 -0.005 
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O.189_ 0005 

fi inn+ - 005 
D.i»u_ 005 

fi 1Q1+0.004 
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O-ly2_ .003 
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O+8'_0.005 
R ls „+0.006 
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D -l°l-0.005 
fi , 7q +0.004 
°-l' y -0.004 
R 17 „+0.006 
O-l'»_0.005 

fi 1 76+ - 006 
D -l' D -0.006 
fi m +0.008 
D ' 1 '°-0.007 
fi 1?4 +0.008 
"■'-' -0.008 
a ,-,,+0.010 
6 - 174 -0.015 
fi 174 +0.004 

D 'i' 4 -o.oo5 

a 1? .+0.003 

D ' 1 '°-0.004 
a m +0.005 
D ' 1 '°-0.004 
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O-l'»_0.007 
a 17 r: + 0.003 
D ' 1 '°-0.004 
fi 17R+0-009 
O-l'»_0.008 

6-i8i±g:8Sg 

6-i86±8;gU 



6.182 1 



f, 1Q9 +0.009 
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b - iyi -0-008 
g 191 +0.007 
o.iyi_ 007 

a mo+°-010 
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fi 1Q4+0.007 

o.iy4_ 00g 
fi 104+0.007 

D - iy4 -0.008 

6.19018:°!; 



6.185: 
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6.185: 



+ 1.5 
1.4 

+ 1.5 
1.2 

1.1 

1.0 
+1.6 

1.4 
+1.6 

1.4 
+1.5 

1.3 
+1.8 

1.5 
+2.0 

1.8 
+2.2 

1.8 
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1.8 
+2.0 

1.8 
+2.0 

1.8 
4 7+ 2 - 2 
4 -'-1.8 
+ 1.7 

1.9 
+2.3 

2.3 
+ 1.1 

2.0 
+ 1.3 

1.3 
+ 1.3 

1.3 
+ 1.3 

1.3 
+77.4 
-2.1 

2.8ii;l 

+2.0 
-1.8 

+2.4 

-2.0 

+2.6 

2.1 



2.6] 
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4.0J 
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3.8" 1 
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4.2" 1 

5.3 4 
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4.5^ 
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2 

2 
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+1.1 



+2.2 

1.8 
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-6.9 
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Q Q+8.8 

y - u -6.6 

= +5.9b 

-4.8 
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R 4^fi+0- 004 
8.43b_ 0004 
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9 7+2.3 
z -' -1.8 
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1 - a -1.0 

9 7+I.6 
z -' -1.4 

2 9+ 17 
z.»_ 1 4 

o 9+I.6 

1.5 
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1.9 

1.6 
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-1.1 

+ 1.8 

1.8 
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1.6 

+ 1.9 

1.7 

1.9 
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1.8 

+2.1 

1.7 

+2.7 

2.0 

+2.3 

1.8 
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-1.4 
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3.8 
4.7 
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4.1 
3.9 
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3.8 
2.7 
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2.2 
1.8 
1.5 
1.6 
1.6 



1.1 

-1.0 

(.0 

1.0 
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0.9 
1.0 
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0.9 



1 O + 0.9 
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1.6+ 09 



-0.8 
+0.9 

0.8 

1 7 +0.9 



1.7 H 



1.9 
3.4 
3.4 



-0.8 

+0.9 

0.8 

1.0 

1.4 
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1.4 



14.0^ 

< 13.2 

13.2i 9 1 3 2 

13.31^ 

.3+8.4 
5 -8.6 
7 Q +4.5b 
'■ y -4.3 
n+5.9b 
8 ' y -6.3 

< 8.7 



10.; 



9.5 



-4.2 b 
-6.4 



< 8.9 

12.8± 7 - 7 

i3.8j:^ 

13.81^ 
1 77 +l7.2 

14 7+14-1 

14t -'-7.9 

H-5±JV 

I0.6j;|;'° b 

< 8.6 
1.0 
1.0 
1.0 
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1.0 

1.0 

1.0 

1.0 

1.0 
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0.008 
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+ 0.008 
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+0.006 
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+0.006 
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-0.008 
+0.009 

0.009 
+ 0.009 
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+ 0.010 
-0.009 
+ 0.010 

0.009 
+0.010 
-0.009 
+0.007 

0.009 
+ 0.008 
-0.009 
+ 0.010 

0.010 
+0.012 

0.013 
+0.019 

0.018 
+0.015 
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+0.020 
-0.019 
+0.019 
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+0.014 
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-2.4 
,+3.3 

-2.7 
- + 1.8 

-1.3 
;+2.4 

-2.3 
3+2.4 

-2.4 
1+3.3 
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7.3J 
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Table 4 

Intcrcombination lines along the phases for 4U1822— 371 a 



phase 



A 



MgXI i 
Flux 



NcIXi 
Flux 



A 



OVII i 

Flux 



1 0.0 

2 0.025 

3 0.05 

4 0.075 

5 0.1 

6 0.125 

7 0.15 

8 0.175 

9 0.2 

10 0.225 

11 0.25 

12 0.275 

13 0.3 

14 0.325 

15 0.35 

16 0.375 

17 0.4 

18 0.425 

19 0.45 

20 0.475 

21 0.5 

22 0.525 

23 0.55 

24 0.575 

25 0.6 

26 0.625 

27 0.65 

28 0.675 

29 0.7 

30 0.725 

31 0.75 

32 0.775 

33 0.8 

34 0.825 

35 0.85 

36 0.875 

37 0.9 

38 0.925 

39 0.95 

40 0.975 
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Table 5 

Fe lines along the phases for 4U1822-371 a 



phase 
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FeXXVI Lya 
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Flux 


a 


A 


Flux 
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1 (upper) 

2 (decreasing) 
3(bottom) 
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5 (rising) 
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i.ao(_ 001 

1 044+0.OO8 
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1 n 36 +0.007 
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a A is in units of 
units of 10" 3 A. 



Flux is in units of 10 °photons s cm 



and 



Table 6 

Ne IX triplet line ratios along the phases for 4U1822— 371 



phase 


A a 
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flux 
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a strongest line in the triplet complex: intercombination line for phase 1,2,3,5; resonance 
line for phase 4. 



Table 7 

Triplet line ratios for two new observations of 4U1822— 371 



Ion 


Hi) 
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(Hr 3 A) ( 
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-1 x10 4 



1x10 4 2x10 
Cycle Number 



Figure 1. Upper left — Light curve for ID 9858 at wavelengths A > 1.5 A; Upper right — The fitted folded phase light curve of 
ID9858: black — observational curve; red — fitted curve. Short black lines and index mark the phase bins for Fe lines analysis. Lower 
left — The 4U1822— 371 partial eclipse timing residuals with respect to the best-fit linear cphcmcris: previous X-ray eclipse measurements 
are at 68% confidence levels — asterisk (Parmar ct al. 2000, and reference therein); previous UV/Optical eclipse measurements at 90% 
confidence — circle; new measurements at 90% confidence( square — ID 671; triangle — ID9076; diamond — ID9858). The best quadratic 
ephemeris from our measurement (dash line) is indicated: dotted line — Parmar ct al. (2000); dash-dotted line — UV/Optical fitting. 
Lower right — Phase folded light curves of ID9858 (black) & ID9076 (red) with example out of eclipse phase bin width shown by blue 
shadows, and the in eclipse phase bin shown by the orange shadows. 
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Figure 2. The phase folded spectra for the five phase bins shown in Figurc[l] 
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Figure 3. Ne IX triplet line spectra in the five phase bins shown in Figure [l] 
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Figure 4. Triplet line spectra Mg XI (left), Ne IX (middle) and O VII (right) for the summed spectra from Obs IDs 9076 and 9858. (See 
Tabled) 
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Figure 5. Line red-shift (left), broadening (middle), and flux (right) as a function of phase for Ly a lines of Si XIV, Mg XII, Ne X, and 
O VIII with 3<7 (solid line) and lc (dotted line) error bars shown. Line broadening is fixed at 0.001A for those points without error bars. 
Green lines mark additional detected line components. 



14 



Ji et al. 



500 - 

o - 

-500 - 
-1000 



MgXI 



0.0 0.2 0.4 0.1 
phase 



0.8 1.0 




'f 


>'f#»: 


."''„' 


MgXI t 



phase 



0.4 O.I 

phase 



1000 



500 - 



-500 

1000 



0.0 0.2 0.4 0.6 
phase 



1500[ 











Ne 


l> 


1000 




l! 1 


% 500 


|.|# 


(fli'llf 

















0.0 0.2 0.4 0.6 O.i 
phase 



7__ 30 

u 25 

en 

,„ 20 



- 


'--/■ 


3!(| T *'< -! 


NelX i i 



0.0 0.2 0.4 0.6 0.! 
phase 



1000 - 

500 



-500 
-1000 



if 



OVII 







ovii i ; 


'e 


T 




OVII 


800 




u 






- 


u 

_ 80 


- 










T 










tr> 






T "' 


600 


T 






III 




- 


£ 60 


j««»'" 


TTt 


T 


4 00 


r >kS' 


""Jl- 


1 


• :; 




* 


- 


a. 40 


- 


,' :l 


3:?§( 


s; ): , T 


?00 


1 


* : 




1 




- 


° 20 


-31 




11 [W< ,. ±i 


'x 






1]- 








,, 


^ 


u: 


,;** »i : 1 





***** 


*** , 


- 


3 o 


,1 


. 


, , 



.2 0.4 0.6 0.8 1.0 
phase 



0.0 0.2 0.4 0.6 0.8 1.0 
phase 



0.0 0.2 0.4 0.6 0.8 
phase 



Figure 6. Line redshift (left), broadening (middle), and flux (right) as a function of phase for the intercombination lines of Mg XI (upper), 
Ne IX (middle), and O VII (lower), with 3o~ (solid line) and la (dotted line) error bars shown, respectively. Line broadening is fixed at 
0.001 A for those points without error bars. Red points mark those phases where the spectra are shown in Figure [71 
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Figure 7. Variations of Ne IX intercombination line at the phases chosen in Figure [B] including the spectra without line detection for 
eclipse (left upper) and phase centered at 0.825 (middle lower). The pink lines mark the zero-velocity line location. The cartoon shows the 
hypothesized structure of the system (based upon the geometrical model of Heinz & Nowak 2001). Each spectrum represents the average 
for the phase bin outlined by the blue shadow on the light curve in the right corner of each figure. 
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Figure 8. Continuum fits to the X-ray spectrum in phase 1. The left panel shows the fit with an exponentially cutoff power law plus a soft 
blackbody and high photo-electric absorption. The power law and black body components are individually shown. The right panel shows 
an alternative solution consisting of a single power law (the cutoff moves beyond the upper energy cutoff of the data) plus highly absorbed, 
partially covered power law with a partial covering fraction of 50%. The overall, ISM contribution to the photo-electric absorption is low. 
The power law and covered power law are individually shown. 
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